Asymptotic performance of optimal state estimation in quantum two level system 
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We derived an asymptotic bound the accuracy of the estimation when we use the quantum corre- 
lation in the measuring apparatus. It is also proved that this bound can be achieved in any model 
in the quantum two-level system. Moreover, we show that this bound of such a model cannot be 
attained by any quantum measurement with no quantum correlation in the measuring apparatus. 
That is, in such a model, the quantum correlation can improve the accuracy of the estimation in an 
asymptotic setting. 
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I. INTRODUCTION 



Estimating unknown quantum state is an important 
task in quantum information. In this paper, we discuss 
this problem by focusing on two typical quantum effects; 
One is the uncertainty caused by the non-commutativity. 
The other is the quantum correlation between particles, 
e.g., quantum interference, quantum entanglement, etc. 
Indeed, the probabilistic property in quantum mechanics 
is caused by the first effect. Hence, it is impossible to 
determine the initial quantum state based only on the 
single measurement. Due to this property, we need some 
statistical processing for identifying the unknown state. 
Needless to say, it is appropriate for effective processing 
to use a measurement drawing much information. There- 
fore, the optimization of measuring process is important 
for this purpose. The second property is also crucial for 
this optimization. This is because it is possible to use the 
quantum correlation between several particles. Hence, we 
compare the optimal performance in presence or absence 
of quantum correlation between several particles in the 
measuring process. This paper treat this comparison in 
the case of two-dimensional case, i.e., the qubit case. 

Estimating unknown quantum state is often formu- 
lated as the identification problem of the unknown state 
which is assumed to be included a certain parametric 
quantum states family. Such a problem is called quantum 
estimation, and has been one of the main issues of quan- 
tum statistical inference. In this case, we often adopt 
mean square error (MSE) as our error criterion. As is 
known in mathematical statistics, the MSE of an estima- 
tor is almost proportional to the inverse of the number 
n of observations. Hence, concerning the estimation of 
quantum state, that of an estimator is almost propor- 
tional to the inverse of the number n of systems prepared 
identically. Especially, it is the central issue to calculate 
the optimal coefficient of it. Moreover, for this purpose, 
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as is discussed by Nagaoka [loj. Hayashi & MatsumotoQ, 
Gill & MassarQ, it is sufficient to minimize the MSE at 
a local setting. (For detail, see section [H]) 

The research of quantum estimation has been initiated 
by Helstromp]]. He generally solved this problem in the 
one-parameter case at the local setting. However, the 
multi-parameter case is more difficult because we need 
to treat the trade-off among the MSEs of the respec- 
tive parameters. That is, we cannot simultaneously re- 
alize the optimal estimation of the respective parame- 
ters. This difficulty is caused by the non-commutativity. 
First, Yuen & Lax Q and Holevo |j| derived the bound 
of the estimation performance in the estimation of quan- 
tum Gaussian family. In order to treat this trade-off, 
they minimized the sum of the weighted sum of the 
MSEs of respective parameters. Especially, Yuen & Lax 
treated the equivalent sum, and Holevo did the generally 
weighted sum. 

After this achievement, Nagaoka 4'j, Fujiwara & 
NagaokajiJ, Hayashi Gill & Massar ^calculated that 
of the estimation in the quantum two level system. They 
also adopt the same criterion. Concerning the pure states 
case, Fujiwara & Nagaoka [Tlj . and Matsumoto[l2j pro- 
ceeded to more precise treatments. 

However, the above papers did not treat the perfor- 
mance bound of estimation with quantum correlation in 
measuring apparatus, which is one of the most impor- 
tant quantum effects. In this paper, we discuss whether 
the quantum correlation can improve its performance. 
For this purpose, we calculate the CR bound, i.e., the 
optimal decreasing coefficient of the sum of MSEs with 
quantum correlation in measuring apparatus, in the sev- 
eral specific model. 

First, as a preparation, we focus on quantum Gaussian 
family, and prove that the above quantum correlation has 
no advantage for estimating the unknown state at sec- 
tion IIIII The reason is roughly given by the following 
two facts. One is the fact that the optimal error with- 
out quantum correlation is given by the right logarithmic 
derivative (RLD) Fisher Information matrix, which is one 
of quantum analogues of Fisher Information matrix. The 
second is the fact that the CR bound can be bounded by 
RLD Fisher Information matrix. 
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Next, we proceed to the quantum two-level system, 
which can be regarded as the quantum analogue of the 
binomial distribution. In this case, as is shown in section 
Ivl quantum quantum correlation can improve the perfor- 
mance of estimation. As the first step, we focus on the 
equivalent sum of the MSEs of respective parameters in 

the parameterization | ( ^ Z x ^ \ w itfi the pa- 

rameter x 2 + y 2 + z 2 < 1. As is discussed in subsection 
IV Al the asymptotically optimal estimator is given as fol- 
lows. 

When the quantum state is parameterized in another 

! ( 1 + r cos 26> re^ sin 29 \ 
way: ^ irk . „„ with the parameter 

2 y re^ 4 ' sm29 l-rcos26 J * 

< r < 1,0 < <p < 2tt,0 < 6> < § , we can divide our 

estimation into two parts. One is the estimation of r, the 

other is that of the angle (9, (f)). 

The estimation of r can be realized by performing the 
projection measurement corresponding to the irreducible 
decomposition of the tensor product representation of 
SU(2), which equals the special case of the measurement 
used in Keyl & Werner[l3j. Hayashi & MatsumotoQ. 
Note that they derived its error with the large deviation 
criterion, but did not treat its MSE. After this measure- 
ment, we perform a covariant measurement for the esti- 
mation of (9, (f>). By calculating the asymptotic behavior 
of the sum of its MSEs of respective parameters, it can 
be checked that it attains its lower bound given by RLD 
Fisher information, asymptotically. That is, this estima- 
tor is shown to be the optimal with the above mentioned 
criterion. Finally, by comparing the optimal coefficient 
without quantum correlation in measuring apparatus, we 
check that using this quantum effect can improve the esti- 
mation error. Furthermore, we treat the CR bound with 
the general weight matrix by a more technical method in 
subsection IV Bl In this discussion, the key point is the 
fact that this model can be asymptotically approximated 
by quantum Gaussian model. 

This paper is organized as follows. First, we discuss 
the lower bounds of asymptotic error in section^ which 
contains reviews of the previous results. In section ITTll 
quantum Gaussian model is discussed. We discuss the 
asymptotic approximation of spin j system by the quan- 
tum Gaussian model in section llVl Using these prelimi- 
naries, we treat quantum two level system in section Ivl 

II. LOWER BOUNDS OF ESTIMATION ERROR 

A. Quasi Quantum CR bound 

Let be an open set in R d , and let S = {pg; 9 £ 6*} 
be a family of density operators on a Hilbert space Ti 
smoothly parameterized by a d-dimensional parameter 
9 = (9 1 , . . . , 9 d ) with the range 0. Such a family is called 
an d-dimensional quantum statistical model. We consider 
the parameter estimation problem for the model S, and, 
for simplicity, assume that any element pg is strictly pos- 



itive. The purpose of the theory is to obtain the best 
estimator and its accuracy. The optimization is done by 
the appropriate choice of the measuring apparatus and 
the function from data to the estimate. 

Let (r(f2) be a a- field in the space f2. Whatever appa- 
ratus is used, the data uo £ f2 lie in a measurable subset 
B 6 tj(ii) of f2 writes 

Pr{^ e B\9} = Pf (B) d = Tr p e M(B), 

when the true value of the parameter is 9. Here, M, 
which is called positive operator-valued measure (POVM, 
in short), is a mapping from subsets B C ft to non- 
negative Hermitian operators in Ti., such that 

M(0) = O, M(O) = / 

M(Q M(B 3 ) (B k n = 0, fc ^ j) 

i=l i=i 

(see p. 53 @ and p. 50 0). Conversely, some appara- 
tus corresponds to any POVM M. Therefore, we refer 
to the measurement which is controlled by the POVM 
M as 'measurement M\ Moreover, for estimating the 
unknown parameter 9, we need an estimating function 9 
mapping the observed data u> to the parameter. Then, a 
pair (9, M) is called an estimator. 

In estimation theory, we often focus on the unbiased- 
ness condition: 

j 9 J M Tr M( dLu)pg = 9 j , V0 e 6. (1) 
Jn 

Differentiating this equation, we obtain 
/ e j {uj)^TtM{duo) Pe = Si (j,fc = 1,2,. 

(2) 

where Sj. is the Kronecker's delta. When (9, M) satisfies 
(JIJ and J5J) at a fixed point 9 £ 0, we say that (9, M) is 
locally unbiased at 9. Obviously, an estimator is unbiased 
if and only if it is locally unbiased at every 9 £ 0. In this 
notation, we often describe the accuracy of the estimation 
at 9 by the MSE matrix: 

V k e ' j (6,M) d = f (9 k -9 k ){fc -9 ] )TrM(dLu)pg. 

or 

tr Vg(9,M)G 

for a given weight matrix, which is a positive-definite 
real symmetric matrix. Indeed, in the quantum setting, 
there is not necessarily minimum MSE matrix, while the 
minimum MSE matrix exists in the classical asymptotic 
setting. Thus, we usually focus on tr Vg(9, M) for a given 
weight matrix. 
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We define classical Fisher information matrix J^ 1 by 
the POVM M as in classical estimation theory: 



•J a ■ — 



wen 



di log — -2-dj log —2- dw 

aw aw 



where di = d/d8 l . Then, Jg is characterized, from 
knowledge of classical statistics, by, 

{Jg'Y 1 = mi{V e {6, M) \ (9, M) is locally unbiased}, 
§ 

(3) 

and the quasi-quantum Cramer-Rao type bound (quasi- 
quantum CR bound) Cg (G) is defined by, 

Cg(G) d = inf{trGV e (0, M) \ 0,M)m locally unbiased}, 
and has other expressions. 
C e {G) 

= inf{tr GVe0,M) \ (9,M) satisfies the condition ©} 

(4) 

(5) 



inf{tr G(J e M )" 1 | M is a POVM on H}. 



As is precisely mentioned latter, the bound Cg(G) 
is uniformally attained by an adaptive measurement, 
asymptotically 0|- Therefore, Cg(G) expresses the 
bound of the accuracy of the estimation without quan- 
tum correlation in measurement apparatus. 



B. Lower bounds of quasi quantum CR bound 

1. SLD bound and RLD bound 

In this subsection, we treat lower bounds of Cg(G). 
Most easy method for deriving lower bound is using 
quantum analogues Fisher Information matrix. However, 
there are two analogues at least, and each of them has 
advantages and disadvantages. Hence, we need to treat 
both. One analogue is symmetric logarithmic derivative 
(SLD) Fisher information matrix Jg-j.k'- 



Jo-j,k d = (Lg-j, Lg ; k}g, 



where 



dpg 



o La- 4 



(X, Y) e = TrpeiX* o Y) = Tr(p e o Y)X* 
XoY d = -(XY + YX), 

and Lgj is called its symmetric logarithmic derivative 
(SLD). Another quantum analogue is the right logarith- 
mic derivative (RLD) Fisher information matrix Jg-j,k : 

Je-j,k = Tr pgLg-k(Lg-j)* = (Lg-j, Lg-k)g 



where 

= peLg-j, (A, B)g d = TrpgBA*, 
and Lgj is called its right logarithmic derivative (RLD). 

Theorem 1 Helstrow^Holevo^ If a vector X = 
[A 1 , . . . , X d ] of Hermite matrixes satisfies the condition: 



(6) 



the matrix Zg(X): 



Z k 6 '\X) d = Trp e A fe A J 
satisfies the inequalities 

Zg{X) > (Jg)- 1 

and 

ZgiX) > (Jg)- 1 . 



(7) 



(8) 



For a proof, see Appendix I A II Moreover, the following 
lemma is known. 

Lemma 1 Holevo^] When we define th vector of Her- 
mitian matrixes Am' 



then 



X 3 m= I V 3 ~(P)M{d6), 



Vg(M)>Zg(X M ). 



(9) 

For a proof, see Appendix IA 21 Combining Theorem 
and Lemma ^ we obtain the following corollary. 

Corollary 1 // an estimator M is locally unbiased at8G 
0, the SLD Cramer-Rao inequality 



Vg(M) > (Jg)- 1 
and the RLD Cramer-Rao inequality 
Vg(M) > (Jg)' 1 



(10) 



(11) 



hold, where, for simplicity, we regard a POVM M with 
the out come in M d as an estimator in the correspondence 
M = Mo§~ 1 . 

Therefore, we can easily obtain the inequality 

trVg(M)G > tv(J e )- l G 

when M is locally unbiased at 9. That is, we obtain the 
SLD bound: 



C$(G) = tr(Jg)- 1 G<Cg(G). 



(12) 



As was shown by Helstrom|2j, the equality l|12H holds 
for one-parameter case. However, we need the following 
lemma for obtaining a bound of Cg(G) from the RLD 
Cramer-Rao inequality. 
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Lemma 2 When a real symmetric matrix V and Her- 
mite matrix W satisfy 



then 



V>W, 



trV > trReW + tr|ImW|, 



where Re W (1m W ) denotes the real part of W ( the 
imaginary part of W ), respectively. 



For a proof, see Appendix IA 31 Since the RLD Cramer- 
Rao inequality (|llfl yields that any locally unbiased esti- 
mator M satisfies 



VGV e (M)VG > VGiJe^VG, 
lemma |21 guarantees that 
tr Vg(M)G 

>trVGRe(J e )- 1 \/G + tr|\/GIm(J ( ,)- 1 \/G|. (13) 
Thus, we obtain the RLD bound: 

Cf(G) d = tr VGRe(Jg)- 1 ^ + tr | VGIm( Jg)" 1 ^] 
< C e (G). (14) 

For characterizing the relation between the RLD bound 
C R {G) and the SLD bound C S (G), we introduce the su- 
peroperator Vg as follows^: 

Pe oV e (X)=i[X,p\. 

This superoperator is called /^-operator, and has the fol- 
lowing relation with the RLD bound. 

Theorem 2 Holevo L r j] When the linear space Tg spanned 
by Le,i, ■ ■ ■ , Le,ci is invariant for the action of the super- 
operator T>g, the inverse of the RLD Fisher information 
matrix is described as 



Jg 1 - Jg 1 + ~Jg 1 DgJ g 1 , 



(15) 



where the antisymmetric matrix Dg is defined by 

De,k,j = (T>e{Lg^),L e . k )g = iTrpg[Le tk ,Le,j}. (16) 
Thus, the RLD bound is calculated as 

C${G) = tr GJg 1 + ^tr\VGJ g 1 DgJ g 1 VG\. (17) 

Therefore, C^(G) > Gf (G), i.e., the RLD bound is better 
than the SLD bound. 

For a proof, see Appendix IA 41 In the following, we call 
the model D-invariant, if the linear space Tg is invariant 
for the action of the superoperator T>g for any parameter 



2. Holevo bound 

Next, we proceed to the non-D-invariant case, in this 
case, Lemma ^ guarantees that any locally unbiased es- 
timator M satisfies 



/ GV e (M)VG > VGZg(X M WG, 
where 

Z k g\X) d =Tv P gX k X3. 
Thus, from Lemma |21 we have 
trVGYg{M)VG 

>Cg(G, X M ) 

d = tr VG Re Zg (X M )VG + tr \VGIm Z e (X M )VG\. 

(18) 

Since Xm satisfies the condition @), the relation (0J 
yields the following theorem. 

Theorem 3 HolevojJ^j: The inequality 



Cf (G) d = min jc e (G,X) 



T ^w XJ = sf ' r " iC:] 



holds. 



Hence, the bound G<f (G) is called the Holevo bound. 
When X satisfies the condition © , the relation 10 yields 
that 

tvGReZg(X) =tiGZ e (X) > tr GJg 1 = Cg(G), 
which implies 

Cf{G) > Cg{G). 
Also, the relation © guarantees that 

VGZ e (X)VG+ VGImZ e (X)VG 
>VGZg(X)VG + VGlmZ e (X)VG > VGJgVG. 
Similarly to (|13|) . the relation © yields 
Cg(G, X) > Cg (G), 

which implies 

Cf{G) > G«(G). (19) 

Moreover, the Holevo bound has another characteriza- 
tion. 

Lemma 3 Let Tg be the linear space spanned by the orbit 
of Tg with respect to the action of T> . Then, the Holevo 
bound can be simplified as 



Gf (G) = min_ \Cg(G,X) 

X:X k £T e I 



Tl w XJ = 6i ■ (20) 
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Moreover, we assume that the D-invariant model con- 
taining the original model has normally orthogonal basis 
(L±,..., L m ) concerning SLD, and the inverse of its RLD 
Fisher information matrix is given by J in this basis. 
Then, the Holevo bound has the following expression. 

Cf(G)= min ftr | VGZj{v)VG\ He{d k \ JW) = 8j\ 

(21) 

where Zj'"* (v) = £ (v k \ J|u J ) and a vector dk is chosen as 



(22) 



Note that the vector v J is a real vector. 



For a proof, see Appendix IA 51 

In the D-invariant case, only the vector L = [Lg d = 
'^2j = i{Jg 1 ) k ' : ' Le-j] satisfies the condition in the right 
hand side (R.H.S.) of (gO), i.e., Cf (G) = Cg(G,L). 
Since Tr pgLgL g = Tr(pgoLg + ^[Lg, pg])Lg, the equation 
(fTo*)) guarantees 

Zg(L) = Jg 1 . (23) 

That is, the equation 

C e (G,L) = Cg R (G) 

holds. Therefore, the equality of (|19() holds. 

Concerning the non-D-invariant model, we have the 
following characterization. 



Theorem 4 Let Si 

dcf 



del' 



C 



{P(e 1 ,...,0 dl ,o,...,o)\(&i, ■ ■ • ,6 

©i} C S 2 =' {pe u ...,e d2 \{0i, ■ ■ ■ ,0 d2 ) C 6 2 } be two mod- 
els such that S 2 is D-invariant. If a vector of Hermi- 
tian matrixes X = [X k ] satisfies the condition 0) and 
X k G< Lg-x, . . . , Lg.d 2 >, then 



C e ,i(G,X) = Cg H 2 (P^GP^ 



(24) 



for any weight matrix G, where the di x d 2 matrix Pjt is 
defined as 



p k dcf Tr dp e xk 



Xd 



(25) 



i.e., is a linear map from a d 2 dimensional space 
to a d\ dimensional space. Furthermore, if the bound 
Cg'2{P x \GPjt) is attained in the model S 2 , the quantity 

Cg,i(G, X) can be attained in the model S 2 - 

Here, we denote the linear space spanned by elements 
v\, . . . ,vi by < vi, . . . , vi >. For a proof, see Appendix 
IA 61 Thus, if the RLD bound can be attained for any 
weight matrix in a larger D-invariant model, the Holevo 
bound can be attained for any weight matrix. 



3. Optimal MSE matrix and Optimal Fisher information 
matrix 

Next, we characterize POVMs attaining the Holevo 
bound. First, we focus on the inequality (|18|l for a strictly 
positive matrix G. if and only if 

V e {M) = B.eZg(X M ) + VG^lVGlmZgiXM^lVG' 1 , 

the equality of (|18|l holds. Thus, the Holevo bound Cg(G) 
is attained for a strictly positive matrix G, if and only if 

V 9 {M)=ReZg(X G ) + VG~ 1 \VGImZ e (X G )VG\VG~ 1 , 

(26) 

where Xq is a vector of Hermitian matrix satisfying 
Cg(G) = Cg(G,Xc)- Therefore, the equation © guar- 
antees that if and only if the Fisher information matrix 
Jf of POVM M equals 



G (VGReZ e {X G )^+\VGlmZg(X G )VG\) VG, 

(27) 

the Holevo bound C e H {G) can be attained by choosing a 
suitable classical data processing. Thus, (|26|l and l|27|) 
can be regarded as the optimal MSE matrix and the op- 
timal Fisher information matrix under the weight matrix 
G, respectively. 

Especially, concerning the D-invariant case, the equa- 
tion (12.'j|) guarantees that the optimal MSE matrix is 

Re(Jg)- 1 +VG~ 1 \VGIm{Jg)- 1 VG\VG~\ 
and the Fisher information matrix is 

VG^VGKe(Jg)- 1 VG+ |\/GIm(J e )- 1 \/G|) 1 VG 

for a given weight matrix G. 

C. Quantum CR bound 

Next, we discuss the asymptotic estimation error of an 
estimator based on collective measurement on n-fold ten- 



def ; 



sor product system 7i®" = Ti <g> • • • ®7i. In this case, we 
treat the estimation problem of the n-fold tensor product 



family S® n d = {pf n = pg ® • • • ® pg \0 e ©}■ Then, we 
discuss the limiting behavior of tr GVg(M n ), where M n 
is an estimator of the family of <S® ra , and Vg(M n ) is its 
MSE matrix. In the asymptotic setting, we focus on the 
asymptotically unbiased conditions l|28|l and l|29|) instead 
of the locally unbiased condition, 



dcf 



6iTrM n (dO)pf n -» & (28) 



<0; k = K, k ( Mn ) = W^ {MTl 



(29) 
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as n — > oo. Thus, we define the quantum Cramer-Rao 
type bound (quantum CR bound) Cg(G) as 



C e {G) 



clef 



mm < lim ntrV e (M")G 



{M n } is asympto- 
tically unbiased 



(30) 

As is independently shown by Hayashi & Matsumoto0 
and Gill & MassarQ, if the state family satisfies some 
regularity conditions, e.g., continuity, boundedness, etc, 
the following two-stage adaptive estimator M™ attains 
the bound C g (G). First, we choose a POVM M such 
that the Fisher information matrix Jg 1 is strictly posi- 
tive for any 9 € 6, and perform it on y/n systems. Then, 
we obtain the MLE 9' for the family of probability distri- 
butions {Pg 1 ^ S 0} based on yjn outcomes ui,..., u/^. 
Next, we choose the measurement Ma, which attains the 
quasi-quantum Cramer- Rao bound Cg(G), and perform 
it on the remaining n — y/n systems. This estimator at- 
tains that C e (G), i.e., trV e (M")G i (%((?). Also, 
it satisfies the conditions (|28(1 and (|29|l . Therefore, we 
obtain 

Cg(G) > Cg(G). 

Moreover, by applying the above statement to the family 
S 0n , we obtain 

nC%{G) > Cg{G), 

where Cg(G) denotes the quasi-quantum Cramer- Rao 
bound of the family S® n . 

In the n-fold tensor product family S® n , the SLD 
Lg, n j and the RLD Lo. n -,j are given as 



where 



Le,n;j — VnL^L Le, n -j — ^/nL^ , 



3 = 1 

X^'fi d = I eg ■ ■ ■ eg) I eg X eg I eg ■ ■ ■ eg J . 
Therefore, the SLD Fisher matrix of S® n is calculated as 

n n 

Tr p^(Le, n ., k o L e , n;j ) = TrpP»(£ L™ o £ L^) 

n n 

1=1 i'=i 

n 

= Y l Trp 9n (L g , k oLg.j)^ 



l-i 

k,j- 



I ■ ■ ■ I ®La._i 



I'-l 



J...J 

n—V 



Similarly, the RLD Fisher matrix of S® n equals the n 
times of Jg. As is shown in Appendix I A 71 a similar rela- 
tion with respect to the Holevo bound holds as follows. 



Lemma 4 Let G^ n (G) be the Holevo bound of 
then 



Cf> n {G) = —Cq(G). 



(31) 



Thus, we can evaluate Cg(G) as follows. It proof will 
be given in Appendix [ 



Theorem 5 The quantum CR bound is evaluated as 

C B (G) > C?(G). (32) 

Therefore, if there exists estimators M n for n-fold ten- 
sor product family such that 

ntrGV e (M n ) -» Cf{G), 

then the relation 

C e (G) = Cf (G) = lim nG e ™ (G) (33) 

holds. Furthermore, if the relation (|33[) holds in a D- 
invariant model, any submodel of it satisfies the relation 

D. General error function 

In the above discussion, we focus only on the trace 
of the product of the MSE matrix and a weight matrix. 
However, in general, we need to take the error function 
g(9, 9) other than the above into consideration. In this 
case, similarly to (|3U|) we can define the asymptotic min- 
imum error Cg{g) as 



Ce(g) 



clef 



mm < lim nH B JM" 

{Af"}~ x [_n->oo 



{M"} is asympto- 
tically unbiased 



where 



clef 



g(9,9)TiM n (d9) 



We assume that when 9 is close to 9, the error function 
g can be approximated by the symmetric matrix G 9 as 
follows: 



9(0,0) 



J2 G U §k 



: )(9 l -9 l ). 



Similarly to subsection lll CI if we choose suitable adap- 
tive estimator M n , the relation K 9 g (M n ) = ±C g {G g ) 
holds. Thus, C e (g) < C (G 9 ). Also, we obtain C e {g) < 
nC%(G 9 ). 
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Conversely, for a fixed 6q, we choose local chart 
at a neighborhood Ug of 8 such that 

g(9 ,e) = ]T Gl^iO) - <P k {9 ))^\0) - <M)), 
fcj 



for V(9 € Ug . By applying the above discussions to the 
family {pe\0 € Ug } 7 we obtain 



C 8 (S) > Cg(G 9 ). 



satisfies the unbiased condition 

Ej(M e ) = 0*. 

Moreover, 1)34(1 guarantees that Tr p g j^M G (d8 1 da 2 ) is 
the normal distribution with the covariance matrix (N + 
\)I + §• Therefore, its error can calculated as follows. 

tr GV e (M 6 ) = tvG((N + hi + §) 



III. QUANTUM GAUSSIAN STATES FAMILY 

Next, we review the estimation of expected parameter 
of the quantum Gaussian state. In this case, Yuen & Lax 
0| derived quasi CR bound for the specific weight matrix 
and Holevo|2j did it for arbitrary weight matrix. This 
model is essential for the asymptotic analysis of quantum 
two-level system. In the boson system, the coherent state 
with complex amplitude a is described by the coherent 

vector | a) = f e~^ L J2^Lo ^fl n )' wnere n) is the n-th 
number vector. The quantum Gaussian state is given as 

dcf 1 f I V I - |Q - C|2 j 
P c ,n = ^ J \a){a\e « da. 



In particular, the relations 



Pom 



-. oo 

+ 1 f< 



N 



N , 

n— 



7V + 1 



|n)(n|, 



hold, where = -^(9 l +9 2 i) and = e J (- e p+e «). 

For the estimation of the family <S^ = {pg^vl^ = 
+ 2 «)}, the following estimator is optimal. Let G 

be the weight matrix, then the matrix G = V det GG _1 
has the determinant 1. We choose the squeezed state 
\(j> G )(<l> G \ such that 



{<P G \Q\<P G ) 1 _ ( 



*g\Q 2 Wg) (4> g \Q o P\^ g ) \ _ g_ 



*e\QoP\4> G ) {<p G \P 2 W G , 



2 ' 



then the relation 



\{<f> G , -h^ 1 + 6 2 i))\ 2 = exp(- k ((G + ir%,jO j 

k,j 



(34) 



holds. The POVM 

MMdd 1 d§ 2 ) Woi,o>\<t> 6 ){<f> d \WZi e 



de 1 de 2 

2tt 



--(N + -) tr G + - tr GVdet GG" 1 
1, 



=(AT+ ~)trG + VdetG. 

For its details, the following theorem holds. 
Theorem 6 Holevo'j] The POVM M G satisfte 



f ${6 l ) 2 M G { d§) J e l 9 2 M G { d§) 

{jo'e 2 M G (de) j(e 2 ) 2 M G (de) 

Q 2 QoP\ VdetG , 



(35) 



(36) 



It is proved in Appendix iBl 

Its optimality is showed as follows. The derivatives can 
be calculated as 



-QQT=-^Pe.N} = WT T(Q-0 ) 



rn i _J_rp 



de 2 



.N 



Therefore, we can calculate as 



N 







</fl 1 Dl)J a 1 — 



1 



(iV + i)- 

-l or 

where we use the relation (|16|) . Thus, since 

^G 



tr 



5^ So 



VdetG, 



the RLD Fisher information matrix is 



iV+i _J/2 
-i/2 7V + i 



Thus, the RLD bound is calculated as 



m F 



G fl fl (G) = (iV + -) trG + VdetG, 
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which equals the right hand side of i|35|l . Thus, from 
l)14p. we obtain the optimality of Mg, i.e., Yuen, Lax, 
and Holevo's result: 

1. 



C e {G) = {N+-)trG + VdetG. 

Furthermore, for the n-fold tensor product model S^™, 
we can define a suitable estimator as follows. First, we 
perform the measurement Mg, on the individual system, 
and obtain n data (9\ , 9f ),..., (9^, 9%). We decide the 

estimate as 9 k d = i?" Ji In this case, the MSE 

matrix equals ±((~N + ~)I + f ). Therefore, Theorem El 
guarantees 



Ce(G) =C e {G) = (N + -) tr G + vdetG, 

which implies that there is no advantage for using the 
quantum correlation in the measurement apparatus in 
the estimation of the expected parameter of quantum 
Gaussian family. 



IV. ASYMPTOTIC BEHAVIOR OF SPIN j 
SYSTEM 



In this section, we discuss how the spin j system 
asymptotically approaches to the quantum Gaussian 
state as j goes to infinity. Accardi and Bachj^, 0] 
focused on the limiting behaviour of the n-tensor prod- 
uct space of spin 1/2, but we focus on that of spin j 
system. Let Jj.i, Jy,2, Jj,a be the standard generators of 
the spin j representation of Lie algebra su(2). That is, 
the representation space Hj is spanned by |j, to),to = 
3, j ~ 1, • • • , -j + 1> -j, satisfying 

J j,3\j,m) = m\j,m). 

dcf 

The matrixes Jj ± = Jj i ± i Jj 2 are represented as 



Jj,+ \j, m ) = Vti - m )ti + m + 1 )\3, m + !> 



Jj-\3,m) = - m + + m)\j, m - 1). 

For any complex z = x + iy, \ z\ < 1, we define the special 
unitary matrix 



U z 



dcf f y/l — \z\ 2 —Z* 



and denote its representation on "Hj by Uj^ z . The spin 



dcf 



coherent vector \j,z) = Uj }Z \j,j) satisfies 
/ / OA 

(j,m\j,z) = 



2j 
j + m 

We also define the state pj yP as 



(i-MT 



Defining the isometry Wj from Ttj to L 2 (R) as 
Wj : \j,m) -> \j - to), 

we can regard the space Tij as a subspace of L 2 (R) . 

Theorem 7 Under the above imbeding, we obtain the 
following limiting characterization 



Pj,v ~* Po,^ 
\z){z\ 



/2j< 



(37) 
(38) 



in the trace norm topology. Moreover, when j goes to 
infinity, the limiting relations 



Trpj, p (a - 
Tr pj, p (a* 
Tr p jtP (Q 
Trp,-, P (P 



l -J ] . + y(a~-±=J J , + )^Q 



/2j 
1 

7W 



Jj,-)*(a* 



1 



/2J 



j x ) 2 -» 



Tr Pj, p Q 2 -> Trp 0j _£_Q 2 
Trpj.pP 2 -> TTp ,_p_P 2 
Tr Pj - p (QoP)->Tr Po ,_P_(QoP) 

Tr Pj , p ((Q- ^J x )oQ)^0 

Vj 

VJ 

Tr Pj> ((P--^J y )oQ)^0 
ZioW, where we abbreviate the isometry Wj 



(39) 
(40) 
(41) 

(42) 

(43) 
(44) 
(45) 

(46) 
(47) 
(48) 
(49) 



V. ESTIMATION IN QUANTUM TWO-LEVEL 
SYSTEM 



Next, we consider the estimation problem of n-fold ten- 

def 

sor product family of the full parameter model 5f u ii = 

{pe =' §i" + Ei=i ^iPH < 1} on thc Hilbert space C 2 , 
where 



1/01 



1 ( -i 



1/1 



471 = 2 I 1 )^ 2= 2 It J' CT3 -2V0 -1 



In this parameterization, the SLDs at the point (0, 0, r) 
can be expressed as 



-k(0,0,r);l — 2<7i, L( ,0,r);2 — 2(72 

1 



^(0,0,r);3 — ^ q -1 



1 — r* 



(2(73 -rJ). 
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Then, the SLD Fisher matrix J(o,o,r) an d RLD Fisher 
matrix J(o,o,r) a t the point (0, 0, r) can be calculated as 



(50) 



We can also check that this model is D-invariant. The 
state pg is described in the notations in Section llVl as 

pe = U e ,4, sin( / >/2 p 1 /2,p(\\8\\)U* i4 , sin4>/2 












i 


1 











n 

!_ r 2 / 






e v sin m. 



where p(r) = j-tl and 

On the other hand, as was proved by NagaokaQ, 
Hayashi^, Gill & Massar[||, in any model of the quan- 
tum two level system, the quasi CR bound can be calcu- 
lated as 



Co(G)= [tiyJJ e *GJ ( 



tr \ VGJrWG 



(51) 



where the second equation follows from the unitary 
equivalence between AA* and A* A. 



A. Covariant approach 

As the first step of this problem, we focus on the risk 
function g covariant for SU(2). Then, the risk function 
R(9, 9) can be expressed by £f(||(9||, 4>), where <j> is the 
angle between 9 and 9, i.e., \(9,9)\ = ||0||||#|| cos cj>. It can 
be divided into two parts: 



g(\\9\\, 11^11^) = /i(l|e||,l|0|l) + /2,||e||.i| fl ||W. 



where 



dcf 



h(\\9\\,\\9\\)=g(\\9\\A\0\\,0) 



{4>) = g{\\ei PU) -sOT.IMI.Q) 



For example, the square of the Bures' distance is de- 
scribed as 

b 2 {pe,P§) = 1 - F(p e ,pg) 



= 2(1- y/^W\l±-\\0\\ 2 -\\0\\\\0\\) 

+ h\&\\ \\e\\{i - cos # 

This risk function can be approximated as 

b 2 (pe,Pe) = \j2 J w( ek - §k W- §l y 

k,l 



Thus, the relations {H}, l|T7|l. and (JSSJ yield that 

C(0Ar)(6 2 ) > C (0,0,r)( b2 ) = — J ' 

Therefore, the covariance guarantees that 
C e {b 2 ) > 



,2s. 3 + 2||0|- 



As another example, we can simplify the square of the 
Euclidean distance \\9 — 9\\ as follows. 



-21 



cos ( 



= (ll<?ll -\\9\\) 2 + 2\\9\\\\9\\(l- cos <j>). 

Concerning this risk function, we obtain 

Ce(I) <Cf{I) = 3 + 2||0|| -||0|| 2 . (52) 

In the following, we construct a suitable estimator for 
the family Sf^ v When we focus on the tensor repre- 
sentation on (C 2 )®" of SU(2), we obtain its irreducible 
decomposition as 

n/2 
j=o or 1/2 

Using this decomposition, we perform the projection 
measurement E n = {E™} on the system (C 2 )®", where 
E'j is the projection to Tij <g> Tt n ,j- Then, we obtain the 
data j and the final state U^^^p^g^U*^^^ <g> 

Pmix,-H„ j with the probability 



Pr,\\e\\(j) 
dcf ( ( n 
1~3 
1- II 



n 

1 + II0II 



7+3 



A - WOLn.-.i + wew.r,. 



n 

1-3 



n 

-3-1 



where PmiyL,n n is the completely mixed state on the 
space Ti. n ,j- Next, we take the partial trace with respect 
to the space Ti. n ,j, and perform the covariant measure- 
ment: 



, sin < 



M'(0, 40 ^ (2j + e** sin ^ sin ||- ^ 

Note that the measure dip dip is the invariant prob- 
ability measure with parameter <f> S [0, 7r),^> S [0,2tt). 



10 



When true parameter is (0, 0, r), the distribution of data 
can be calculated as 



For a proof of (J5SJ, see Appendix IE II 

Therefore, the average error can be approximate as 



Tr ft)P M^,V) 



1 - p 2 ^ 



T l-(l-p)sin" 



2 / 4tt 



where p = p(r). 

Finally, based on the data j and (0, t/>) , we decide the 
estimate as 



-i 2j ,j 2 2 7 -o 2 7 

? = — cos sin 0, 8 = — simp sin (j>, 8 — 



COS( 



Hence, our measurement can be described by the POVM 

M" v d = {M j (<j), if)) <g> J Wn ^ } with the outcome (j, <j>, tp) . 

Next, we discuss the average error of the square of the 
Euclidean distance ||t9 — 8\\ 2 except for the origin (0,0,0). 
For the symmetry, we can assume that the true parame- 
ter is (0, 0, r). In this case, the average error of \\9 — 9\\ 2 
equals 

£ P„ ) ((|- r ) 2 + 2,-!^), (53, 
j=o or 1/2 \ x ' I 

where 



F, 



clef 



27T f 7T 



o Jo 



(1 - cos <t>) Tr p. i-r M 3 (0, ip) d<t> dip 

Jl i_|_ r 



(l-cos^) ( 1- (l-p)sin^- 



2j 



sin i 



2(l-p2j+i) < /_ 1 v- ^ 2 ' 2 

_2 (1 + (2j + l)p 2 J+ 2 - (2j + 2)p 2 J +1 ) 
(2j + 2)(l-p)(l-p 2 J+i) 

Thus, for any fixed p < 1, we have 
1 1 



2j 



1 - p j + 1 



o(i)p 



2j 



(54) 



Using the above relation, the first and second terms of 
(|53|l can be calculated as 



\8-8\\ 2 T r M? ov (d8)pf n 



9N 1 ,2 1 - r) 4,1 

= 3 + 2r - r 2 - - (-i ^ + W. 

v y n r 1 + r ?i 2 

Combining this and l|52|) . we obtain 

C e (/) = Cf(/)=3 + 2||0||-||0|| 2 . 



(56) 



The average of the square of the Bures' distance also can 
be calculated by the use of the relations l|54|) and l|55|) as 



Up, 1 - VT^Jl - (^) 2 - 2 ±r + — ^ 
2 \ V n n n J ' !+>■ 

-c 1 

2 V 

1 



'- (1 - r 2 ) + - (1 + r) — 

4(l-r 2 ) 1 i+oV^')- 



where we use the following approximation 



1- V / T T ^ 2 A/l-(-) 2 --r = 



1 (2j 



n n 2(1 — r 2 ) \ n 
for the case when -2 is close to r. Thus, 

As a byproduct, we see that 
2j 

► r as n — > oo 

?i 

in probability -P n , r . 

Next, we proceed to the asymptotic behavior at the 
origin (0, 0, 0), In this case the data j obeys the distribu- 
tion P n , : 

^*M(r-iK-"i-.)) w+i >- 



As is proved in Appendix IE 21 the average error of the 
square of the Euclidean distance can be approximated as 



^,M[ 2 {) (57, 
\ n J n y/irny/n n z 



n/2 

, s 2j 2 Ar 2 x . 1 . 



,j=o or 1/2 

z — ' \ n 

^=0 or 1/2 x 

= (1 - r 2 )- - 2 _^A\ + 0((1 _ r 2 ) n/2 ) _ 

n r n z 



(55) 



Since 

9 - (0, 0, 0)|| 2 Tr M™ v (^)pf ^ )0) = ]T P n ,o(j) 

(58) 

we obtain (7(0,0,0) CO = ^(o o o)(0 = ^> *- e -' * ne e Q ua tion 
JSS) holds at the origin (0, 0,0). 
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On the other hand, by using l|51(l . the quasi quantum 
CR bound can be calculated 



C e {I) = (2 + Vl-ll^ll 2 ) 2 = 5 - ||0|| 2 + 4 V /T^pjp. 

(59) 



(3 + 2||0| 



2 ) = 2(1- 



Since 5 - 

||0||) + 4-^/1 — ||#|| 2 is strictly greater than in the mixed 
states case, using quantum correlation in the measuring 
apparatus can improve the estimation error. 

Remark 1 The equation 154|) gives the asymptotic be- 
lt can 



havior of the error of M J ((f), tp): Fj_ p = nzz^n ■ 
be checked from another viewpoint. First, we focus on 
another parameterization: 



M j (z)dz A = (2j 



The equation 138fl of Theorem guarantees that the 
POVM M^{^=) goes to the POVM \z)(z\. Thus, the 

equation 1371) guarantees that its error goes to with the 
rate nj^j ■ This fact indicates the importance of ap- 
proximation mentioned by Theorem Indeed, it plays 
an important role for the general weight matrix case. 

Remark 2 One may think that the right hand side 

(R.H.S.) of (|55|l is strange because it is better than 

(1 — r 2 ) — , i.e., the error of the efficient estimator the 

. n . . . 

binomial distribution. That is, when data k obeys n- 

binomial distribution with parameter ( i=£ , i±£ ) and we 
choose the estimator of 9 as k/n (it is called the efficient 
estimator), the error equals (1 — r 2 )-^, which is larger 
than the right hand side of (|55|) . However, in mathe- 
matical statistics, it is known that we can improve the 
efficient estimator except for one point in the asymptotic 
second order. In our estimator, the right hand side of i|55|) 
at r = is given in Q57J1 ■ and is larger than (1 — r 2 ) — . 



B. General weight matrix 

Next, we proceed to the general weight matrix. For the 
SU(2) symmetry, we can focus only on the point (0, 0, r) 
without of loss of generality. Concerning the RLD bound, 
we obtain the following lemma. 

Lemma 5 For the weight matrix G = ^ ^ ^ ^ , the 
RLD bound at (0, 0, r) can be calculated as 

C (o,o,r)( G ) = trG- r 2 .s + 2rVdetG, (60) 

where G is a 2 x 2 symmetric matrix and g is a 2- 
dimensional vector. 



For a proof, see Appendix IE 31 The main purpose of this 
subsection is the following theorem 



Theorem 8 Assume the same assumption as Lemma\^ 
then 



C (0fi ,r)(G)=Cf 00r) (G)=trG-r 2 s + 2rVdetG. 

(61) 

Furthermore, as is shown in Appendix lE 41 the inequality 

C e (G) = C?(G) <C e (G) (62) 

holds. Thus, using quantum correlation in measuring 
apparatus can improve estimation error in the asymptotic 
setting. 

As the first step of our proof of Theorem [SJ we 
characterize the MSE matrix attaining the RLD bound 
C£ i0jP) (G). The matrix 



satisfies Vq T > J^ QQ * and 



dot f 1 + rVdetd ■ G~ x 



1 - r 1 



tr GV e r = tr G - r 2 s + r tr VdetG • G _1 G 



= tr G - r 2 s + 2rVdetG = Cf Q r) {G). 

Thus, when there exists a locally unbiased estimator with 
the covariance matrix V^, the RLD bound G^ 0r - ) (G) 
can be attained. 

In the following, we construct an estimator M n locally 
unbiased at (0, 0, ro) for the n-fold tensor product family 
S f ®n such that riV (ofi , ro )(M n ) ->V G . In the family Sg™, 
the SLDs can be expressed as 

j 

for k = 1, 2, and 

First, we perform the projection- valued measurement 
E n = {Ej 1 }. Based only on this data j, we decide the 
estimate of the third parameter 8% as 



Jn,r 



dr 



(63) 



where 



ri log P„, r (j) 
dr 
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Then, we can easily check that this estimator 0% satisfies 
the following conditions: 



Tr 



dpf n 



d8 k 



1 k = 3 



(E^) = l0 fc = 1)2 

(0,0, r) j ^ 



(64) 
(65) 



for Z = 1, 2. The definition of Of(j) guarantees that 



The definition guarantees the equation (|65l) and the equa- 
tion l|64|) for k — 3. The rest case can be checked as 
follows. The derivative of pg with respect to the first or 
second parameter at the point (0, 0, r) can be replaced 
by the derivative of U x +i V p(ofi yr )U* +iy with respect to x 
or y. Since the probability TrU^pf^Uf^yM, is 
independent of x + iy, we have 



dTrpf n E] 



for k = 1,2, 



(66) 



which implies Q64JI in the case of k = 1,2. 

Next, we take the partial trace with respect to H nj -, 
and perform the POVM V*M ( j( dx 1 dx 2 )Vj on the space 
TLj. After this measurement, we decide the estimate of 
the parameters Q 1 , O 2 as 



3' r \ X 



where 

dcf f Ti(pj, P o 2J jA )VfQV j Tr(p 3 , p o 2J 3 , 2 )V*QV J 
- { I '''/'..;■ ° 2J i , 1 )^/Py j Tr( Pj , p o 2J j , 2 )V*PV j 

As is shown in Appendix IE 51 the relations 



Tr 



(0,0, r) 

<96> fc 



(67) 



^^0?0,r)(©(^ M i,G(^)® J «n.i) -0 (68) 



hold for I — 1,2, k — 1,2,3. Therefore, we see that our 
estimator (O 1 ,^ 2 ,^) is locally unbiased at (0,0, r). 
Next, we prove that its covariance matrix V n satisfies 



v , I + rVfdG-i V =y 1. (6g) 
1 — r } n • n 



Using the equation i|E9|) in Appendix IE 51 we have 



T ^foA,) (© -r)M Ae (d»)) ®/«„, 



As is shown in Appendix IE 61 the above value can be 
approximated by 



Jn.l = (1 " ^) 



1 1 - r l 1 



(70) 



In order to discuss other components of covariance ma- 
trix, we define the 2x2 matrix V • 



j,G,r- 



[Vjy = [Tr PjtP J x k x%M d {dx)V*]. 
By use of Theorem El this matrix can be calculated as 



Tr p,.,,\,Q 2 \: Tr p jtP Vj (Q o P)V/ \ VdetG ^ 
Tr ft , p ^(QoP)V* Tr ftiP ^P 2 U/ J + 2 ^ ' 



then the covariance matrix of the estimator q on the 
state Pj,p( r ) is 

B i,r Y j,G,r( B hr) T - 

Theorem \7\ and l|3l)l) guarantee that 



as j —>■ oo. Hence, 



Thus, the covariance matrix of our estimator (f? 1 ,^ 2 ) 
equals 

£ P n Aj)B-^A B lr f = (! + rV / det(5 • G- 1 )- 



because the random variable — converges to r in proba- 
bility P ntr . Thus, we obtain JdS}. 

Concerning the origin (0, 0, 0), we can prove 



-J M - -J 

„ J (0,0,0) ^ J ' 



(71) 



which will be proved in Appendix IE 71 Therefore the 
RLD bound at the origin (0, 0, 0) can be attained. Then, 
the proof of Theorem |S| is completed. 
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Holevo bound in submodel 



VI. DISCUSSION 



Since the Holevo bound is attained in the asymptotic 
sense in the full model, Theorem [5] guarantees that the 
Holevo bound can be attained in the asymptotic sense in 
any submodel S — {pe^lv £ 6 C R d }. In the following, 
we calculate the Holevo bound in this case. Since the 
Holevo bound equals the SLD Fisher information in the 
one-dimensional case, we treat the two-dimensional case 
in the following. First, we suppose that the true state is 
P(o,o,r)- Without loss of generality, by choosing a suitable 
coordinate, we can assume that the derivatives can be 
expressed as 



del' 



dp8( n ) 



def 9p9(r)) 



9=(0,0,r) 



cos 0(72 + sin 



r 2 cr 3 , 



0=(O,O,r) 

J. In the above assumption, we have the 



drj 2 

where < <\> < 
following theorem. 

Theorem 9 Assume that the weight matrix G is param- 
eterized as 



G 



91 .92 

92 53 



31 



< the Holevo bound G* 0r) (G) of 



When 

VdctG 

the above subfamily can be calculated as 



C(0,0,r) (^) 



tr G + 2r cos ^Vdct G - r 2 sin 2 . 



and can be attained only by the following covariant matrix 



V G = / + r cos 0VdetG • G -1 



r 2 sin 2 



(72) 



Otherwise, the Holevo bound C% r )(G) and the covari- 
ant matrix Vg can be calculated by 

detG 



trG 



.91 



V G = J 



cos" 



sm 



/ COS 


*)' 


\ sin 


t>) 


3 2 2 

— 92 


_ 92 

91 
1 


91 



(73) 
(74) 



For a proof, see Appendix IE! 

On the other hand, the equation (|51|l guarantees that 



G (0 ,o,r) (G) = (tr VG) 2 = tr G + 2Vdet G 

in this parameterization because Jg — I. Since we can 
verify the inequality tr G + 2\/det G > G^ r ^ (G) in the 
above two cases, we can check effectiveness of quantum 
correlation in the measuring apparatus in this case. 

The set {V G |detG =1} represents the optimal MSE 
matrixes. Its diagonal subset equals 

j . ( rt~~ x cos d> — r2 a i r > 2 









rt cos < 



< t < 



cos< 



r sm <p 
(75) 



We proved that the estimation error is evaluated by 
the Holevo bound in the asymptotic setting for estima- 
tors with quantum correlation in the measuring appara- 
tus as well as for that without quantum correlation. We 
construct an estimator attaining the Holevo bound. In 
the covariant case, such an estimator is constructed as 
a covariant estimator. But, in the other case, it is con- 
structed based on the approximation of the spin j system 
with sufficient large j to quantum Gaussian states family. 

It is also checked based on the previous results that 
the Holevo bound cannot be attained by the individual 
measurement in the quantum two-level system. That is, 
using quantum correlation in the measuring apparatus 
can improve the estimation error in the asymptotic set- 
ting in the quantum two-level system. 

Since the full parameter model of the quantum two- 
level system is D-invariant, its Holevo bound equals the 
the RLD bound. Thus, its calculation is not so diffi- 
cult. However, a submodel is not necessarily D-invariant. 
Hence, the calculation of its Holevo bound is not trivial. 
By comparing the previous result, we point out that this 
model is different from pure states model even in the 
limiting case r — > 1. 
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APPENDIX A: PROOFS OF THEOREMS AND 
LEMMAS IN SECTION El 

1. Proof of Theorem 111 

For any complex valued vector b = [bj] and we define 
a complex valued vector a — [a-*] — Jg l b and matrixes 
= f J^j X^bj and = f J2j Le-jaK Since the assump- 
tion guarantees that {X^, L£) = (b, a) , Schwarz inequality 
yields that 

(b\Z e (X)\b)(b\J^\b) = (b\Z e (X)\b)(a\J e \a) 
={X^X^){L B ,L S ) > | (b, a) | 2 = \{b\J^\b)\ 2 . 

Therefore, we obtain 

(b\Z e (X)\b) > (blJ-'lb), 

which implies Q. Similarly we can prove JSJ). ■ 
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2. Proof of Lemma [D 

For any complex valued vector b = [bj], we define ma- 
trix X b,Ai = Ej X M h r Sincc 

§,b)M(d§)=X SM , 

we obtain 

(b\V e (M)b) - (b\Z e (X M )\b) 

(§,b)*{e,b)M(de)-xi M x lM 

((9,b)-X ?M rM(d§)((9,b)-X VM )>0, 
which implies JJJJ • B 

3. Proof of Lemma [2] 

def 

Since the real symmetric matrix T = V — Re W sat- 
isfies 



T > Im Wi 



we obtain 



trT > min{trT'|T' : real symmetric, T > ImW} 
= tr|ImW|. 

Therefore, 

tr V > tr Re 14 7 + tr T > tr Re VT + tr | Im W\ . 



Next, we define a linear map £ from C d to as follows, 

i 

then its inverse £ _1 and its adjoint £* are expressed as 
d 

C- 1 :X^J2(Jg 1 )^(L 9 . !k ,X) e 

k=l 

C* :X»(L e .j,X) e . 
Thus, the map Jg can be described by 

£* o (/ + i-X^)" 1 o £ = £* o P Te (7 + i-Vg^Pr, o £, 

where Pr„ is the projection to Tg. Since is invariant 
for X> e , 

(P Te (7 + l ip e )- 1 P T J- 1 = P Te (7 + i\v e )P Te . 
Therefore, the inverse of Jg equals 

C- 1 o (P Te (7 + l l -Vg)- 1 P Te y 1 o (£*)-* 

=£- 1 oP Te (7 + l il? e )Pr e o (r- 1 )*, 
which implies 

(^- 1 ) M = E( J fl" 1 ) M ^'( 7+ 4^ )L ^' )e(Je ~ 1) '' ,J - 



5. Proof of Lemma l3l 



4. Proof of Theorem El 



Sincc 



6 r 

P6» o Lg.j = pgLg-j = pgo Lg-j + ~[Lg-j,pg\ 



-P8 ° ( Lg.j + -Vg(Lg.j) ) , 



we have 



[I + i-Vg){Lg,j) = Lg.j, 

which implies Lg.j = (7 + i^T>g)~ 1 Lg-j. Since 
Iff (^)* = (Wpa. ^ have 

Je-k,j = Tr PeLg-k{Lg.j)* = Tr(Lg-j)* pgLg ;k 

= Tr ^-Lg, k = Tr(p e o Lg.j)Lg. k = (Lg.j, Lg ; k)$ 

= (Le, J ,(I + t^Vg)- 1 Lg, k )g. 



Let P be the projection to Tg with respect to the inner 
product ( , )g, and P c be the the projection to its orthog- 
onal space with respect to the inner product When X sat- 
isfies the condition ©, {P(X k ),Lj) e = (X k ,P(L J ))g = 
(X k ,Lj)g = S k . Thus, P{X) = [P{X l )\ satisfies the con- 
dition ©. Moreover, 

Tr pgP(X k )P c (X j ) 

= Tr fpg o P{X k ) + ±\p, P(X k )]\ P C (XJ) 



= Tr f pg o P(X k ) + -poV e (P(X k ))j P c {Xi) 
= Tr {pg o (p{X k ) + ±Z> fl (P(X*)))) P c (^) 

= (P(X k ) + l -Vg(P{X k )),P C {Xl))g = 0. 

Thus, we obtain 

Zg(X) = Zg{P{X)) + Zg{P c (X)) > Zg(P(X)), 
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which implies that 

VGReZg(X)VG+\VGImZ e (X)VG\ 

>Vgz 9 (x)Vg > Vgz 8 (p(x))Vg. 



7. Proof of Lemma l4l 

Let Tg be the linear space spanned by the orbit of the 
SLD tangent space of S® n . Since any element X of Tg 
satisfies 



Since the matrix VGIm Zg(X)y'G is imaginary Her- 
mite matrix, \>/GIm Zg(X)y/G\ is real symmetric ma- 
trix. Therefore, Lemma [21 guarantees that 

Zg(X) > Zg(P(X)), 

which implies i|2U[l . 

Next, we proceed to a proof of H21|) . Since the basis 
(Li, . . . , L m ) is normally orthogonal concerning SLD, the 
equation l|15|) guarantees that 



(n) 



Tr pgL k Lj = Tr p e L k o Lj + - Tr pg[L k , Lj] = S k j - i-D, 



.1 

'2 
(Al) 



e,k 



Hence, when we choose the vector v k — (v k , . . . v^J sat- 
isfying that X k = J^j v j L ji 



Tv^X k = -Re(d k \J\v^) 
TrpX k X j = (v k \J\v j ). 
Therefore, we obtain (|21|l . 

6. Proof of Theorem |H 



(A2) 
(A3) 



There exists d\ x d,2 matrix O such that ^JP^GP-g 
OVGP^ and T = I dl . Since X k = J2fli P \, X L V 



Ce,i(G, X) 

-- tr VG Re Zg {X)VG + tr \y/G Im Zg {X)VG\ 
-- tr VGPji Re Zg {L) T P^ VG 
+ tr \VGPf lmZ e {L) T P^^/G\ 



= tr OVGP^ Re Zg{L) T P^VG O 

+ tr \0VGPg Im Z g (L) T P^VG 7 0\ 
=C e , 2 (P^GP^L) = C% 2 (P%GP jt ). 

Let {M™} be a sequence of locally unbiased estimators 
of S 2 such that Tr Vg (M n ) PZ GPg -» Cl 2 { T P^GP^). 

Next, we define an estimator M n> = f (Pg,M n ) on 5i 
satisfying locally unbiasedness condition at 9. Its covari- 
ance matrix is V$(M n ') = V r e (M") J P|. Hence, 



trVe(M"')G = trVg{M n )P^GP^ 
^{ T P i GP jt ) = C e ,i{G,X). 



V^iVg o • • • o £><,(X)J = V^Z> 9 o • • • o P«(lW), 

the equals 

{V^x(")|XeT e }. 

Furthermore, the vector y/nX^ = [y/n(X')^] satisfies 
Cg(G, y/nX^) = nCg(G, X). 

Therefpre, Lemma [3] guarantees that 
,j = ^6 ■ 
Cf' n {G) 

= min_ \C (G, y/fticW) (VnL$, ^(X ] ) {n) )g = 5{ } 
= min_ {nC^G.X) n{Lg, k ,X^)g = S 3 k } 

X:Xi£T e < > 

where we put Y = —X. Therefore, we obtain (|31|l . ■ 



: min_ <^ -Cg(G,Y) 

Y:Yi£Tg l n 



8. Proof of Theorem |5] 

Lemma n guarantees that 

Yg{M n ) > Zg{X Mn ). 

Since the vector Y M = (3* d = ^.^(M)" 1 )]!; of 
Hermitian matrixes satisfies 

{\fnLg. % ,Yl {n )g = 5?, 

^(Im») - A 9 (M")Z e (y M «)^(M"), 

the relations 

tr GWg{M n ) 

>tr VG-ReAg(M n )Zg (Y M « )Aj(M n )VG 

+ tr \VGImAg(M n )Zg{Y M ™)Aj{M n )VG\ 

>C^ n (Aj(M n )GAg(M n )) 



1 



=-Cf (M™)G?j4e(M™)) 
n 



hold. Taking the limit, we obtain 
* lim C 

n — >oo 
:Cf(G) 



lim ntrGV e (M")> lim Cf {Aj (M n )GAg(M n j) 

l— >oo n — >oo 



which implies l|32|) . 
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APPENDIX B: PROOF OF THEOREM [6] 

Let E be the joint measurement oiP®I + I®P and 
Q®I-I ®Q on the space L 2 (R) ® L 2 (R) . As was proved 
in Holevo |]f , the POVM Mg satisfies 

TrM d (dxdy)p = TrE(dxdy)(p®\(l> d ){<t> <s ,\). (Bl) 

Thus, 

Tr x 2 Mg,(dxdy)p = Trx 2 E(dxdy)(p ® |</> G )(</> G |) 
= Tr(Q 2 ®/ + /®Q 2 )(p® |0 G )(0 G |) = (TrQ 2 p) + G 1 ' 1 , 

which implies equation i|36|) regarding the (1,1) element. 

Since (</> G |P|0 G ) = (</> G |Q|</> G ) = 0, Concerning (1,2) 
element, we have 

Tr xyMg(dx dy)p — Tr xyE(dx dy){p ® |</> G )(0 G |) 
= Tr(Q oP®I + I®QoP-P®Q + Q®P) 

■{p®\<t> & ){4>G\) 
=(Tr(Q o P)p) + (0 G |Q o P|0 G ) = Tr(Q o P)p + G 1 ' 2 . 

We can similarly prove equation (|36J) for other elements. 

APPENDIX C: PROOF OF THEOREM 

First, we prove JHZJl. Since 

Pj-p ~ Po,y?^ 

2j+l 2 J 00 

= 1 l p 2 J - + i ( 1 -p)Ep"l n )H-(i-p) E p»h 

^ n=0 n=2j + l 



we have 



Ilft-.P-Po.^ll 

p 2j + l ^ 



^rrd-^E^ + d-p) E 



n=0 



n=2j + l 



<- 



1 - p 2 ^ 1 



-> 0, 



which implies (|37(l . Next, we prepare a lemma for our 
proof of (jg) . 

Lemma 6 Assume that a sequence of normalized vector 
o" = {a™}°^ and a normalized vector a = {ai}°l Q satis- 
fies 



Proof: For any real number e > 0, there exists an 
integers Ni such that 

oo 

£>| 2 < C . 

i=JVi 

Furthermore, we can choose another integer such that 

Ni-l ATi-l 

io?-a<i a <c, e ikt-n 2 i < e > V ™>^ 

j=0 i=0 

Hence, we have 

oo Ni — 1 /iVi-1 \ 

^K| 2 = 1-^K| 2 <1- ^k| 2 -e <2e. 

»=M i=0 V i=0 / 



Therefore, 



JVi-i 



E i< - «<i 2 < E 1°? - a *i 2 + 2 E (Ki 2 + n 2 ) 



<e + 2(2e + e) = 7e. 

Then, our proof is completed. 
We can calculate \j, -|=) as 



E 



2j _ 



2j-re 

|2\ -V- 



2j 



Its coefficient converges as 



2j 

2j-n; VV27 



1 - 



(2j) 



(2j-n)!(2j)' 

,„,2 a ™ 



2j 



2.7 

2\ -™/ 2 / 



27-5 

|2\ -V 



1 



2j 



as j — > oo. 



Thus, Lemma |3] guarantees that 



\z)-\j,-=) 

v 2 j 



which implies l|38(l . 



then 



a, — > a,- as n — > oo, 



Ek-«.i 2 -°- 

2 = 



J A+ |n) - v^V2j-" + l|" - 1) (n = 1, • • ■ ,2j) 
J i)+ |0)=0 



Jj,_|n) = v 7 ™ + 1 V 2 J - "I" + 1) (n = 0, . . . ,2j - 1) 
J,-,_|2i)=0 
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I Tr Pj, P Q 2 - Tr p p Q 2 1 + \Trp j<p P 2 -Tr Po » P 2 



L — p 

~i2 



i-p 
1 - P 2 n 



2j 

tE 



' 2j - n + 1 



2j 



l))V|n-l)(n-l| 



| Tr( Pj , p - p 0)T ^)Ql + I Tr(p JiP - Po.^P'l 

p 2j+l 



Since the inequality 1— yl — x < ^fx holds for < x < 1, 
we have 



2j 



< 



oo 

-(1-p) ^ p»(n|)(Q 2 + P 2 ) 

n=2j'+l 



n=0 



- 1 _ p 2j+l 27 P 

^ n=l J 



1 - p 2 J+ 

Tr(l-p) ^ p»(n|(Q 2 + P 2 ) 



<a-f)E^ 



2 I p pi 1 • /" 



1 2j 2j (l-p)3 



0. 



n=2j + l 

: I £^ fT (l_p)^(2n+l) 



n=0 



which implies 1|39[) 
1 



(a* 



2j J j,-)PjA a * - ^Wj J - ! - r 



(1-p) ^ P n (2n+1) 

n=2j'+l 

2p 



2j-l 



n=0 



+ T 2 ^T^ + 1 )V^|2 J + 1>(2, + 1|. 



i-p (i-pr 

oo 

+ (l-p) E P"( 2 «+!) 

n=2j + l 

->0 as j — ► oo, 



Since the inequality 1 — — x < y/x holds for < x < 1, 
we have 

Tr 0* - ^Jj,-)PjA a * ~ ^J J ^-)* 



because Y^=iP n {^ n + 1) < °°- Thus, we obtain ll4"3|) 
and iJHJ. 

|Trp J)P (QoP)-Trp 0)rf _(goP)| 

,2j+l 2 J 



< 



< 1-P (n + l)n „ 

-!_pM+i 2^ 2 j P 

1 n=0 J 



1-P 



l - p 2 J 

oo 

(l-p) J] p»<n|)(QoP) 



<a-p)E 



2j (l-p) 3 l-p 2 J 
which implies l|40|l . Since 
1 



" r (2i + l)V J '^0 



n=0 J ^ 

1 — P 1+p 1— P 



n=2j+l 
2j+l 2 J 



Since 



(g--^j K ) 2 + (p--^j y ) 2 

=(<* - -i? J i,+)*( a - 7^ J J>+) 
+ (« - 4+)( a - -i J j,+)*> 



the relations l|39|) and H4U|) guarantee the relation 1)41(1. 
Also, we obtain ijl^jl. 



^1 Tr i! 2j+1 (i-p)E>>)M(Q ° m 

oo 

+ |Tr(l-p) 53 p n \n)(n\(QoP)\. 

n=2j + l 



\{Q 2 + P 2 )<QoP< l -{Q 2 + P 2 ) 



|Tr(l-p) 53 p n \n){n\(QoP)\ 

n=2j+l 

<Tr(l-p) f] p"|n)(n|i(g 2 +P 2 )^0 

n=2j+l 
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Similarly, we ca show 



Proof: The Fisher information Jg equals 



n2j + l 



2j 



y ^ — n 



Thus, we obtain 14511 . By using Schwarz inequality of the 
inner product (X, Y) i— ► Trpj jP (X oY), we obtain 



< Tr Pj , p (Q - -^J x ) 2 Tt Pj , p Q 2 . 

Vj 

Thus, the relations l(41j) and 143(1 guarantee the relation 
igSJ. Similarly, we obtain 1071) - (gHJ). 



dlogp 9 (^i)pg(^ 2 |cji) 



dlogpe(^i) . dlogp e (o;2|wi) 



d0 



<20 

dlogp e (o;i) 



d0 



+ ^ pe(w2|wi) 



d(9 



d0 



d0 



APPENDIX D: USEFUL FORMULA FOR FISHER 
INFORMATION 



In this section, we explain a useful formula for Fisher 
information, which are applied to our proof of <|7(J[) and 
CU- Let S = {p e (u}i,uj2)\0 G 6 c M} be a family 
of probability distributions on Sli x f2 2 . We define the 
marginal distribution pe(wi) and conditional distribution 
as 



However, the second term is vanished as 

dlogp e (w 2 |wi) 



dO 



d\ogp g {uJi) dlogp e (wi) ^ dp e (w 2 |wi) 

2^ Pfl(wi) ^ — 2^ 



de 



w 2 GS2 2 



=0. 



Pe(wi) = ^ Pe{u x ,u)2), pe{u2\ui) 



dcf w 2 ) Thus, we obtain (|D1|) . Moreover, we can easily check 



Then, the following theorem holds for the family of dis- 
tributions S, S\ = f {pe(uJi)\9 G 6 C M}, and S Ul = f 
{p e (w 2 | Wl )|e G 6 CM}. 

Theorem 10 The Fisher information Jg of the family S 
satisfies 



Je = Ji.e + ^ Po{ui)Jui, 



(Dl) 



where J\fi is the Fisher information of Si and J ull fi is the 
Fisher information of S u , 1 . Moreover, the information 
less has another form: 



Ju lt e = E w(w 2 |o;i) 
W2sn 2 



dlogp d (uJ2,U!i) 



de 



E Pe(^2\uJi) 
\w 2 e^2 



d\og P g(uJ2,Ul) 

d6 



(D2) 



Thus, the average Ylu efi can be regarded 

as information loss by losing the data u!%. 



APPENDIX E: PROOFS FOR SECTION [V] 



1. proof of (1551) 



The L.H.S. of l|55|l can be calculated as 



i/2 



2j 



j=o or 1/2 

fc=0 ^ 

/ [ n /2] /p *. 

=4 g(r) 2 + V P n , r (n/2 - A) ^ - 2g(r)- 



where q(r) d = ^p- Since the probability 
Pn,r{% — k) has another expression: P n ,r(f — k) = 



K^-^-i))^ 1 -^))"^ 1 ! 1 -^) 



n—k 
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we can calculate the expectations of k and k 2 as follows. Furthermore, every term appearing in the above equation 

is calculated as 



[n/2] or [n/2]-i 



k=0 



k( n )q(r) k (l-q(r)y 



-k 



k=0 



)q(r) k (l-q(r)y 



+ 0((l-r 2 )™/ 2 ) 
--np + 0{{1 - r 2 ) n l 2 ) 

[n/2] Or [n/2]-l 



]T . Wr)*(l-g(r))' 

: ^fc(fc-l)f^ 9 (r) fc (l- g (r))"- fc 



-k 



k=0 v 7 

+ 0((l-r 2 )™ /2 ) 
=n{n-l)p 2 + 0{(\-r 2 ) n ' 2 ). 

[n,2] 

^ ■fri.rf'g — *0 Note that (1 — r 2 ) > Using there formulas, we 

k=0 obtain 



V2] 

g(r) fe (l-g(r))"- fe+1 



fe=0 





.G 








— r 




+ r 



n/2N 



[ ™ /2] 1,2 jU 

g(r) 2 + £ P n , r (n/2 - - 2 g (r)- 

* — ' n 

k=0 

1 - r 2 1 1 - r 1 



- , , o +0((l-r 2 )"/ 2 ), 

[n/2] . . 4 n 2r n 2 



I^(fc(fc-l) + fc)(")g(r) fc (l-g(r)r 

, [n/2]-l , . 

- J2 (^-l)+3HCL(rf +1 (l- 9 ( r )r 



which implies Ij55(l . 



+ (;((! — ") ) 2 ' P roofof <E3 



[n/2] 



The left hand side of l|57[) is calculated as 



fc=0 



£ *; (C) -(*!,)) «rtO - «r>r-» £ ? (d - j) - (f -" - 1)) <2J + »(» 



1 [n/2] /„\ 1 AA/ s /n-2-0\ 



2 



r ^ \k U=o 

/2X / 



O 



1-r 

1 + rJ I klo 



, (n-2fc-l)(n-2fc-2) 2 (El) 
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When n is even, (n - 2(f) + l)(n - 2(f)) 2 = . Then, 



the above value are calculated 
-l 



n 2 2 



'.On 



k=0 



(n-2fc + l)(n-2fc) 2 



- ('')(^-2fc- l)(n-2fc-2) 2 
fc=0 "' 



When n is odd, (n - 2[f ] - l)(n - 2[f ] - 2) 2 = 
Then, the above value are calculated 

\fc=0 v 7 

-EffcV- 2 ^ 1 )^-^- 2 ) 2 ) 

fc=0 ^ 7 / 



n 2 2 r 



E 



(n-2fc) 2 -8(n-2/s) + 4 



The first term is calculated as 
(n-2fc) 2 



- y 

2" ^ U 



k=0 v 7 

= - ■ 6 • 4n 2 • — = 3n. 



2 An 

-1 /'n-lA 



k=0 



n—l tn-l\ 
k ) 



Since £L=o (V) = Em« ( 
1 



t j6(,-2-) 



we have 



n x I n—l 



2 n+l 



E 

fc=0 
n—l 

E 



fc=0 



m [u on+l 



2 n+l 
1 

2 



fc=0 

n\ 1 



A- 

71 — 1 



n 
= 4 

n\ 1 



— 7 2 n+1 



2 n+l 



' ' E (?)6(n-2fc) 2 -8(n-2fc)+4 ] . (E2) 



n 2 2" \^ \ k 

\k=0 



The first term is calculated as 



2™ ^ \ k 

k=0 



2 n+1 ^ \kj 2 n 

1 „ 1 
=- • 6 • 4ri 2 = 3n. 

2 4n 



Since EltV 1 (V) = XX[Wi (V)- we have 



2™ ^ U / ' = 2" ^ I k 

fc=0 v 7 fc=0 v 

n_1 / i\ / 1 

// / // - - I \ ;/ / // — I 



k=0 v 7 v 12 



n n (n — l 

ra 

- 2 



4 2«+! V [f] 



and 



2"^Ui 2"+ lZ ^U/ 2 

fc=o v 7 fc=o v 7 



1 l*J 
2 n ' 

k=Q 



Thus, 



n / On+l 
2/ Z 



(n - 2k) + 4) 



1 fn\ 

2" E L J (~ 

fc=0 v 7 



fe=0 



2 ™+i 



,n / n 



1 n 1 

™y2"+ T_ 'ir ^2 



" / On+l < 
2/ Z 



=(- 8 -- 4 )(a)2^T) + 2 - 



Since (")^ 



2 -^, we have 



Since (p])^r — \ <|E2|) can be approximated as 



R,H.S.oME21 



1 



1 / n n f n — l 



_^3n + (-8n + 4)- + 16^-^1 



If n n - 1 



2 3 4V2 2 



n n^/7r[f ] n 2 n y/ifny/n n 2 



3. Proof of Lemma [5l 

First, we parameterize the square root of G as 



i_4(2l + J-)f n 
n V n n 2 ; U 7 2 n + l 



1.2 3 
) + — = - 



4V2 1 



/7T ri\ n r 



VG= i T 



^4 a 
a T t 
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where A is a 2 x 2 symmetric matrix and a is a 2- 
dimcnsional vector. 



^(0,0,r) (^) 



trG-r 2 s 2 +rtr 



A a 
a T 



-i 

1 




A a 
a T 



By putting J : 
term as: 



-1 

1 



, we can calculated the second 



tr 



= tr 



A a 
a T 



-i 

1 




A a 
a T 



(detA)J AJa 
(AJa) T 



=2^(det A) 2 + \\AJa\\ 2 = 2^/{detA 2 ) + (Ja\A 2 \Ja) 
=2^det(A 2 + \a){a\), 



where the final equation can be checked by choosing a 
HI ^ r: - i- 

obtain (fHTT|> - 



basis such that a — 



Since G = A 2 + \a)(a\, we 



4. Proof of (|62|> 

First, we focus the following expressions of Cg(G) and 
C e {G) 

C e {G) = (tr y/y/GJf l VGj (E3) 

C?(G) = tr VGJ^VG + tr ftVGJ^DoJ^VG]. 

(E4) 

When we put the real symmetric matrix A = 



VGJg 1 v / G and the real anti-symmetric matrix _B = f 
2VGJ f r 1 D g J- 1 VG, the relation 



A + iB > 0. 
Here, we dragonize A as 

' a 
1 | & 

0c 



(E5) 



(E6) 



with a > b > c and c > 0, where the final strict inequality 
follows from G > 0. Since \B\ is a constat times of a two- 
dimensional projection P. Hence, 

PAP + iB = P(A + iB)P > 0. (E7) 

If we regard PAP as a two-dimensional matrix, tr \B\ < 
2v det PAP. Thus, by considering the maximum case of 
the minimum eigen value of PAP, we have 



tr\B\ < 2\fri> 



(E8) 



Therefore, 

Cg(G) - Cf(G) = (tr VA) 2 - (tr A + tr \B\) 
>2 ( Vab + Vbc + ^fca - Vabj = 2 (Vbc + ^/caj > 0. 



Since 



we have 



5. Proofs of (Hill! and (HHJ) 



VJU.r'KXVj \ Q p k k Z\ 



Tr{p j}P o 2 J A fc ) / ^M i)S ((i0) = 4 



for fc, Z = 1, 2, where 



dcf 



M J - ) g( d6) = V/Mg o Bj\d0)V v 

Since the matrixes o * 2 (2.7^3 — rl) and pj p are 
diagonal and all diagonal elements of V*QVj and V*PVj 
are 0, we have 

Tr (^p ° 7^72 ( 2J ^3 " rI))V;QYi = Tr p jiP v;QVi = 
Tr( Pj> o I ^ 3 (2J Jt3 - rI))V;PV 3 = Tr p^V* PV 3 = 0. 
Thus, 

Tr(p„ P ° 7^-3(2^,3 - rl)) f 6 l M 3 6 { d§) 
= Tr Pj> / § l M j6 (d§)=0 (E9) 
for Z = 1,2. Therefore, 

= E (?) Tr (<°^ ° 2 / ^ g( ^) = 4 

for fc, i = 1,2. For the k — 3 case, the above quantity 
equals 



1 — r Jr 



= 0. 



Furthermore, we have 



= E ^VO') Tr ^.p / ® lM iM d§ ) = ( E1 °) 



for 1 = 1,2. 
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6. Proof of (JToJ) Thus, 



First, we focus on the following equation 
n 



1-r 2 



*o?o.r) f©i^(2Ji.a-ri)®J«. rf j 



2 

(2m-r) 2 . 



Then, applying Theorem 1 11)1 we can see that the differ- 
ence " 2 — J n ,r equals information loss. Thus, 



Jn,r ^ ^ Pn.r (j) ,r ; 



1 _ r 

where 



^ dcf ^ 1 -p(r) -_ m / 2m_r x 2 

J ' r ~ ^ l-p(r) 2 J+ lW I 1-r 2 



)2 
T - T 2 " ~ J "' r ' r 2 (l-r 2 )' 

y j-yw P{ rr m ( 2{m ~ j) ) 2 

V 1- ^ r) DfrV- m 2(m ~ j) > i 

i_«( r )2j+i^ ^ 1-r 2 
m=-j ) 

4(1 (p(r)+p(r) 2 



m=-j 



(1 - p(rfi+i)(l - r 2 ) 2 \ (l-p(r)) 3 
_ (l+p(r))p(r) 2 J +1 

(4j + (2j) 2 (l-p(r))Mr) 2 ^+ 1 \ 

(i-pW) 2 ; 

4(1 -p(r)) 2 



(1 -p(r) 2 J+ 1 ) 2 (l -r 2 ) 2 

'p(r)(l -p(r) 2j ') 2jp(r) 2 J' +lx2 



(l-p(0) 2 1-P(r) 
+ 0(p(r) 2 ^ 



r 2 (l-r 2 ) ' " which implies (25 . 
I 



7. Proof of (f7l|> 

For the covariance of the POVM M^, v , the Fisher information matrix J, Q is a scalar times of the identical 
matrix. We apply Theorem 1521 to the family of probability distributions p r (j, 4>, ip) = f TrpS, n r ^M^{(j), ip) ® ^-h„ >3 = 
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Pn,r(j) Tr Pj^ p M 3 ((f>, ip). Then, we calculate the Fisher information: 

On the other hand, Applying Theorem 1 101 to P n ,r(J)(j-> m \Pj,p\3> m )> we have 



n 



1 _ r 2 



j J m=-j 



m) 



Thus, 



J ^= T ^r 2 -Y. P ^)[J [ d ^ ^ ^ ) 2 Tr Pj , p Mi {ff,, 



- 2^ ( J r J 0,m\p jt p\j,m)\. 

m=-j ^ ' / 

In the case of r = 0, Since Tr pj >p M^ (<j>, i/,) = (2j + 1) ^ffi+i (i 1 ^)'' we obtain 



o Jo 
Jo 



± ) Tl " P]-P(r) M J (0, V) # dip 

* ( d^j^r(^) 23 ) 2 . !_ p(r) f 1 + rcM ^ 



dr I (2j + 1) l-p(r)^i I 1 + r J 



o Jo 



Its first and second terms are calculated as 

2^ 2 



Jo 



2 " f f rfl °g( r+r ) I , , l-p(r) fl + rcos^V 3 sin0 fl, n2 

' ( 2 J' + 1 ) 1 ^I 2+1 , , „ ~7Z i ^diP= -(2jcos^) 2 sin 



dr | w ; 1 - p(r) 2 J' +1 V 1 + r / 4tt 7 2 

2 J 



G?r Jo Jo ^ r 1 — p(r) 2 J+ 1 \ 1 + r / 47r 

* 1 

— 2j cos <p sin </> <i</> = 0. 



Since we have 



J 

Therefore, if the relation 



\ - / d(j,m\ Pj , p \j,m) \ 2 . 4 ..4. 4\/2 ^2 

^ I Jr ) ^MPi,p\3,m) = l3{3 + ^), P rh0 (j)-j =^=V^ + ^ (Ell) 



m=-j 
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holds, we obtain 17111 . Hence, in the following, we will 
prove i|Kll|l . 



2 -y-(( 



i-i 



2j(2j + 1) 



4gM(*H-i> l( "- 2fc+1)( "- 2t) 



2 1 /n 



3 2" V 



(n-2-0 + l)(n-2-0) 



2 ' z3 (?)(n-2fc + l)(n-2A) 



3 2™ \ ^— ' \ fc 

\fc=0 



[f]-i 



fc=0 



53 ( )(n-2fc- l)(n-2fc-2) 



(E12) 



When n is even, (n - 2(f) + l)(n - 2(f)) 2 = . Then, 
the above value are calculated 



53 P ».o(i) 
j 

4rfE (?Vfi-2fe+l)(n-2fc) 



fc=0 



53 (")(n-2Jfc-l)(n-2A!-2) 



When n is odd, (n - 2[f ] - l)(n - 2[f ] - 2) 2 = . Then, 
the above value are calculated 

53 p «>oO')oj 



(E(3(n-2* + l)(n-2*) 
\fe=0 v ' 



k=Q 



53 (")(n-2Ar-l)(n-2jfe-2) 



4\/2 



30? V 3' 



which implies ilEllfl . 

8. Proof of Theorem I 



We focus on the full parameter model with the deriva- 
tives at the point (r, 0, 0) 

d P dp dp 2 

oei = ai ^ = a2 '^ = {1 ~ r )as (E13) 



96 )2 " Q 3 



as a D-invariant model. In this case, the SLD Fisher 
information matrix is the identity matrix. Thus, we can 
apply (|2*TJ) of Lemma Hence, by putting 








ch = cos < 



we obtain 



C^(G) = min Ut\VGZj(v)VG\ Re(d k \J\v j ) 
where 




dcf 



Hence, from the condition 

(d 3 \J\v k )=S k j . 

Then, v 1 and v are parameterized as 

v 1 =L\ — t sin 4>Li + t cos 4>L 3 

v 2 =(— s sin <f> + cos 0)^2 + (s cos + sin (f))!/^. 

The matrix Zj{v) can be calculated as 

1 + t 2 ts — ir(—s sin <f> + cos <^>) 

is + ir(—s sin + cos 0) 1 + s 2 

Thus, the quantity tr |v / G-^j(w)v / G| equals 



tr G + 5i (i + —s) H s 2 + 2r| cos0 - sin<£s|Vdet G. 

(E14) 



< 



In the following, we treat the case of 

The minimum value of (|E14|) equals tr G+2r cos <p det G- 



r sin (f>gi which is attained by the parameters t = 
— 33-s s = rgl =£ . Thus, the discussion in subsubsec- 

tion III B~3l guarantees that the Holevo bound is attained 
only by the following covariance matrix 

ReZe(X) + VG~ 1 \VGImZ e (X)VG\VG~ 1 



1 + t 2 
ts 1 



ts 



l (f) + cos ^iVdet GG- 1 
(E15) 



=R.H.S. of O- 

In the opposite case, the minimum value of (|E14|1 equals 
R.H.S. of l|73(l . which is attained by the parameters 
t = — — s.s — cos f . Substituting these parameters into 

gi ' sin or- 

(|E15)l . we obtain O. ■ 
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